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ABSTRACT 

We investigate the topology of the new Point Source Catalogue Redshift Survey 
(PSCz) of IRAS galaxies by means of the genus statistic. The survey maps the local 
Universe with approximately 15000 galaxies over 84.1 per cent of the sky and provides 
an unprecedented number of resolution elements for the topological analysis. For com- 
parison with the PSCz data we also examine the genus of large N-body simulations 
of four variants of the cold dark matter cosmogony. The simulations are part of the 
Virgo project to simulate the formation of structure in the Universe. We assume that 
the statistical properties of the galaxy distribution can be identified with those of the 
dark matter particles in the simulations. We extend the standard genus analysis by ex- 
amining the influence of sampling noise on the genus curve and introducing a statistic 
able to quantify the amount of phase correlation present in the density field, the am- 
plitude drop of the genus compared to a Gaussian field with identical power spectrum. 
The results for PSCz are consistent with the hypothesis of random phase initial con- 
ditions. In particular, no strong phase correlation is detected on scales ranging from 
10 /i~^Mpc to 32 h~^Mpc, whereas there is a positive detection of phase correlation at 
smaller scales. Among the simulations, phase correlations are detected in all models at 
small scales, albeit with different strengths. When scaled to a common normalization, 
the amplitude drop primarily depends on the shape of the power spectrum. We find 
that the constant bias standard CDM model can be ruled out at high significance 
because the shape of its power spectrum is not consistent with PSCz. The other CDM 
models with more large-scale power all fit the PSCz data almost equally well, with a 
slight preference for a high density rCDM model. 

Key words: galaxies: clusters: general - cosmology: observations - cosmology: large- 
scale structure of Universe 



1 INTRODUCTION 

All-sky redshift surveys of galaxies selected from the IRAS 
catalogues have had a dramatic impact on our understand- 
ing of Ihe large scale struct ure and dynamics ot the local 
Univeri 



our Local Group motion or the local velocity field as deter- 
mined from peculiar velocity studies, with the velocity field 
inferred from the IRAS galaxy density, h as provided deter- 
minations o f 0,'''^ /b on very large scales ( Rowan- Robinson 



et al. 1990|; |Kaiser et al. 199l|; Pekel et al. 1993| ; [Nusser & 



e. inet^jjuu, survey (^^wrcnccetal. l^^o[) provided pavis 1994| ; ^trauss et al. 1992bD . These dynamical stud 



early evidence that the large-scale clustermg ot galaxies was 



incompatible with the standard cold dark m atter model ( Ef 



3tathiou et al. 199C; 



Saunders et al. 1991). Comparison of 



ies were in general compatible with Q, = 1 cosmogonies with 
only small deviations of the bias parameter from 6=1. This 
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suggested that IRAS galaxies are reasonably faithful trac- 
ers of the mass field, as compared to e.g. optically selected 
galaxies. The main advantage of IRAS galaxy redshift sur- 
veys is, however, the uniform all-sky coverage as a result of 
the IRAS survey planning and the reduced Galactic obscu- 
ration in the far infrared. IRAS galaxy surveys were also 
able to probe cosmologically significant volumes efficiently 
due to their broad selection function. 

The two most significant IRAS redshift surveys to date 
employed two different strategies. The QDOT survey con- 
sists of 2387 galaxies randomly selected at a rate of 1 in 6 
to the full depth of the IRAS Point Source Catalog (PSC), 
0.6 Jy, thus exploring the maximum vol ume at the expens e 
of sampling density. The 1.2-Jy survey (Fisher et al. 1995 ), 
an extension of the 2-Jy survey (Strauss et al. 1992a), com- 
prises 5321 galaxies fully sampled from the PSC but to a 
shallower depth. 

The sparse sampling strategy of QDOT was ideal for 
low order statistics such as the galaxy power spectrum 
whereas the 1.2-Jy survey was more useful in measure- 
ments of higher order statistics wh i ch require higher sam 



1995; E 



plin g densities ( Fisher et al. 1994; [Musser, Dekel fc Yahi] 



ouchet et al. 1993). Some important studies such 



as the inferred local velocity field, the convergence of the 
dipole and the topology greatly benefit from both high sam- 
pling densities and cosmologically significant volumes. For 



thes e re asons we have constructed the PSCz ( Saunders et al _ 
1994 |), ^ complete redshift catalogue of 15000 galaxies to the 
fuU depth of the PSC (0.6 Jy). 

The construction of the parent sample for the PSCz 
was similar to that of the Q MW IRAS Galaxy Catalogue 
( Rowan- Robinson et al. 1990| ) with substantial efforts made 
to improve the sky coverage, uniformit y and completeness. 
Details are provided in Saunders et al. (1997). Optical coun- 



terparts were identified using digitized plate material from 
APM and COSMOS. Galactic sources were excluded on the 
basis of IRAS colours and literature searches. Existing sur- 
veys and literature material provided around two thirds of 
the redshifts, allowing the full redshift catalogue to be con- 
structed with realistic amounts of telescope time. We have 
now measured around 5000 new redshifts to provide a fi- 
nal catalogue with redshift completeness of 98 per cent; the 
observing pro gramm e and data reduction are described in 
Keeble et al. jl997| ). 

There has been extensive discussion in the literature 
on the ability of topological studies to discriminate between 
different models of structure formation; not only between 
inflationary models (with random-phase initial conditions) 
and non-Gaussian models of structure formation, like cosmic 
strings or textures, bu t also between models with difi'erent 
dark matt er content (Kerscher et al. 1997 ; Matsubara & 
Suto idgi; [CoUey, Gott fc Park 199e|). From these previous 



studies, it is clear that the genus statistic is a robust dis- 
criminator as long as low-noise data are used. Our use of 
very large simulations of CDM variants in conjunction with 
the PSCz survey is ideal to reveal the power of this statis- 
tic in discriminating between different models of structure 
formation. 

Topological studies have also been aimed at drawing 
conclusions on the shape of the power spectrum. Moore et al. 
(1992) concluded that the amplitude of the genus curves on 
large scales was inconsistent with the predictions of a con- 



stant bias standard CDM model and that the power spec- 
trum was best fit w ith a power law with index n — — 1. 
Vogeley et al. (1994) showed that the values for the am- 
plitude of the genus curves were again inconsistent with a 
constant bias standard CDM model but consistent with an 



open CDM model, and Protogeros & Weinberg (1997) con- 



cluded that the topology of the 1.2-Jy redshift survey was 
best fit by a power spectrum with index n — —1. However, 
the quantitative significance of these results was not up to 
what one would wish, given the limitations of the data they 
used. 

This paper forms part of a sequence describing the first 
results of the PSCz survey: Saunders et al. (1997) describe 
the cosmography and an alvse the counts-in-cells statistics; 
Rowan- Robinson et al. (1997) examine the converge nce o f 
the density and velocity dipole; Sutherlan d et a l. (1997) 
discuss the power spectrum; Tadros et al. (1997) perform 
a spherical harmon ic an alysis of redshift space distortions, 
while Keeble et al. (1997) compare the radial and transverse 
components of the correlation function. Here we analyse the 
topology of the PSCz density field: Section H discusses the 
topological statistics applied while Section describes the 
construction of the PSCz density map. We present results 
for PSCz in Section ^ and give an account of the theoretical 
models used to compare with the data in Section H Finally, 
we present our conclusions in Section |^. 



2 TOPOLOGICAL METHODS 
2.1 Genus statistics 

The statistical tools currently used in order to confront 
theories of structure formation with observational data go 
well beyond the simple two-point correlation function. The 
search for such 'higher-order' information is based on the re- 
alization that the two-point correlation function (or equiva- 
lently the power spectrum) exhausts the statistical content 
of a system only when this has a Gaussian nature. It is 
clear, however, that non-linear gravitational clustering will 
inevitably introduce non-Gaussian features, even allowing 
for random-phase initial conditions. 

In order to meas ure d epartures from Gaussianity, Gott, 
Melott & Dickinson (1986) suggested studying the topology 



of isodensity surfaces and quantifying it via the genus. Sub 
sequently the genu s statistic has b e en applied to a number 
of redshift surveys (|Gott et al. 198£|; [Moore et al. 1992| [Park 



Gott fc da Costa 1992^ [Vogeley et al. 1994|; [Rhoads, Gott'fc 



subject to several t heoretical investigations ( riamilton, Gott 



fc Weinberg 19861; IWeinberg, Gott fc Melott 1987 



Weinberg fc Gott 1988|; [Park fc Gott 199l|; [Beaky, Schcrrer 



Postman 1994; Protogeros fc Weinberg 1997 ) and has been 



Melott 



fc Villumsen 1992[ iGott, Gen fc Ostriker 1996|; [Matsubara & 



Yokoyama 1994; Matsubar^'ir^okoyama 1996^ ^a^subara 



199(: ; Matsubara fc Suto 1996). The genus has also been con- 



sidered in the more general framework of Minkowski func- 
tionals ( Mecke, Buchert fc Wagner 1994 ) . 

Given an isodensity surface S, the genus G of that sur- 
face can be defined as 

G{S) = # of holes — # of isolated regions -I- 1, (1) 
i.e. a spherical surface has genus 0; a torus has genus 1, 
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whereas a distribution of A*' disjoint spherical surfaces gives 
rise to G = -(TV - 1). 

A more formal definition of the genus can be given by 
means of the Gauss-Bonnet theorem, which relates the cur- 
vature of the surface to the genus. Let k = (aia2)^^ be the 
local Gaussian curvature of a two dimensional surface S, i.e. 
ai and 02 are the two principal radii of curvature. Then, the 
Gauss-Bonnet theorem states that 

G(S) = — - / kdA + 1. (2) 
Js 

Defining S as the boundary surface between two re- 
gions above and below a threshold density pt we can calcu- 
late G{S) and plot it as a function of that threshold, thus 
obtaining the genus curve of a given density field. 

What do we expect the genus curve to look like? If a 
high threshold is selected, only a few very dense and isolated 
regions will be above this density value and the genus is 
negative. If a low threshold density is chosen, only a few 
isolated voids are identified and, again, the genus is negative. 
On the other hand, for a threshold around the mean density 
value, one expects, in general, that the isodensity surfaces 
will have a multiply connected structure that resembles a 
sponge, with a resulting positive genus. 



The genus curve defined as a function of u has the 
remarkable property that it is invariant during the linear 
growth of density fiuctuations in the Universe. The density 
contrast Sp/p grows and contours of p change, but contours 
of 1/ are fixed. This is a very important property as by mea- 
suring the genus curve on scales where non-linear growth 
has yet to occur we can recover information about the pri- 
mordial density field and, eventually, be able to distinguish 
between Gaussian and non-Gaussian initial conditions. 

As an alternative to labeling the isodensity surfaces 
with the number of standard deviations from the mean, it is 
common to parameterize them with the fraction /vol of the 
survey volume above the given density threshold or, equiv- 
alently, by the number Vf given implicitly by 



/vol = 



(27r) 



1/2 



-t^/2 



dt. 



(6) 



Note that the two parameterizations coincide for the 
case of a Gaussian random field. The volume fraction param- 
eterization has the disadvantage that it is less sensitive to 
skewness in the density probability distribution. In addition, 
it does not discriminate between Gaussian distributions and 
any other distribution which is a one-to-one transformation 
of the former, for example, a lognormal distribution (Coles 



dom field. In this case, the statistical properties of the den- 
sity field are completely described by the power spectrum 
and so the genus will not provide any additional information. 
Moreover, underdense or overdense regions are statistically 
indistinguishable. This means the genus curve is symmet- 
ric about the mean, which is characteristic of the so-called 
'sponge-like' topology. Parameterizing the threshold density 
by the number v = [p — Pmcan)/o"p of standard deviations 
from the mean density, and introducing the genus per unit 
volume g( S) = (G(S) — ^)/V , one finds for a Gaussian ran- 
dom field ([Hamilton, Gott fc Weinberg 1986|): 



Th ^ simplest case occurs, of course, for a Gaussian ran- fc Jones 1991 ). This is because the contours of constant vol 



2x -v^ /2 



Here the amplitude 



N 



3/2 



(2vr) 



is a constant which depends via its second moment 
/ PP(fc) d^fc 



{k^) = 



J P{k)d^k 



(3) 



(4) 



(5) 



ume fraction are invariant under such mappings, as well as 
the topology of the iso densi ty contours themselves. 

As Vogeley et al. ( 1994 ) argue, this invariance property 
is indeed desired, if one is not interested in the positive skew- 
ness developing during gravitational collapse, but rather in 
the topological properties of the initial density field. By us- 
ing the volume fraction to label isodensity contours, the 
genus curve becomes independent of the one-point proba- 
bility distribution function (PDF). In this way the topologi- 
cal analysis is not mixed up with the one-point PDF, which 
can be better studied by other means. Subscribing to this 
philosophy we will hereafter use the volume fraction param- 
eterization. 



2.2 Calculation of the genus curve 

We want to utilize the genus to study the topology of the 
observed galaxy distribution, which comes in the form of a 
point set. As a first step we therefore need to adopt a method 
to compute s uitabl e surfaces from the point distribution. 
Mecke et al. (1994) assigned to each galaxy site a ball of 



on the power spectrum P{k) of the (smoothed) density field. 

Hence, the genus curve of a Gaussian random field ex- 
hibits a universal w-shape. Only the genus amplitude de- 
pends on the shape of the power spectrum via equation (|^, 
but it is independent of the normalization of P{k). 

For a non-Gaussian density field, different topologies 
are expected. For example, the genus curve can be shifted 
towards a 'meatball' topology or towards a 'Swiss-cheese' 
topology. In a meatball topology, the genus curve peaks at 
a negative value of u because isolated structures dominate 
over an extended range at positive densitiies, whereas in a 
Swiss-cheese topology, the genus curve peaks at a positive 
value of since the topology is dominated by empty voids 
up to a larger density value. 



radius r and examined the genus of the union set of these 
spheres. We will employ the more widely used approach of 
constructing smoothed maps of the PSCz density field (see 
Section 3) and considering isodensity contours for the genus 
statistic. 

Our code to calculate the genus curve of a given 
smoothed density field is ba sed on the algorithm first pro- 
posed by Gott et al. (|l986|) and on the CONTOUR code 



by Weinberg (1988). We have written a new implementa- 



tion in C, able to compute high resolution genus curves very 
quickly. Here a major increase in speed was achieved by sort- 
ing the discrete density field first in order to be able to find 
the threshold values pt for the required values of /vol at no 
computational expense. 

Given a particular surface of constant density, we ap- 
proximate it by a network of polygonal faces. This procedure 
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does not change the global topology of the surface as long 
as the grid size is much smaller than the smoothing length. 

When such a polygon network is used to a pproximate a 
compa ct surface o f genus G, it can be shown (pott, Melott 



Dick 



nson 1986) that the genus is given by 



^A = 47r(l-G), 



(7) 



where Di = 360° — ^ . Vi is the angle deficit at each vertex 
and Vi are the angles around the vertex. For example, in 
a cube, which approximates a compact surface of genus 0, 
there are three squares around each vertex and, thus, Di — 
360° — 3 X 90° = 90° at each of the eight vertices, giving 
. Di — 720° = 47r, as expected. The reason for this is that 
the curvature is compressed into ^-functions at the vertices. 
Parallel transport arguments show that the integral of the 
5-functions of Gaussian curvature over the infinitesimal area 
of a vertex is just equal to the angle deficit at that vertex. 

In practice we divide space into cubic lattices for sam- 
pling either galaxy or simulation data; we then calculate the 
genus of a given isodensity surface by approximating it by a 
network of square faces, and by adding up the angle deficits 
of all vertices. 



2.3 Genus related statistics 

In order to measure departures of the observed genus curve 
from the random-phase shape we use t he genus meta- 



statistics introduced by Vogeley et al. ( 1994 ) . They consist of 
an appropriately defined amplitude, width, and shift of the 
genus curve. Additionally we consider the amplitude drop of 
the genus curve compared to the random phase expectation. 



2. 3. 1 Amplitude 

We measure the amplitude of the genus curve as the ampli- 
tude A'^ of the best fit random phase genus curve by mini- 
mizing -x^ in the range — 1 < i/ < 1. If the underlying density 
field is sufficiently close to Gaussian the amplitude pro- 
vides direct information on the shape of the power spectrum. 
It should be noted, however, that the amplitude is system- 
atically biased high by shot noise due to finite sampling of 
the density field, as will be discussed below. 



2.3.2 Width 

The second of the meta-statistics is the width of the genus 
curve, which refers to the range of v over which the genus is 
positive, i.e. 



Wv = v+ — V- 



(8) 



where 1/+ and V- are the first zero crossings of the genus 
curve, right and left of the origin. The width Wv measures 
whether the examined density field is more or less sponge- 
like than a random-phase field, which has Wv = 2. In order 
to improve the reliability of the determination of v+ and 
V- we first reduce the noise in the genus curve by boxcar 
smoothing it with a filter of total width Av = 0.2. This 
smoothing is only applied for this statistic. 



2.3.3 Shift 

In order to distinguish between a Swiss-cheese and a meat- 
ball topology, we use a measure for the shift of the peak of 
the genus curve. This shift may be quantified via 



(9) 



where g{v)ohB is the observed genus curve an d q(i/)fit is 



the c orresponding best-fit random-phase curve (Park et al 
1992). We evaluate Av in the range 



1 < < 1. A value of 
Al/ > indicates a bubble-shift, whereas Av < means a 
meatball-shift. 



2.3.4 Amplitude drop 

With the above meta-statistics, any positive detection of a 
deviation of the genus curve from the random phase shape 
signals the presence of non-Gaussian features in the density 
field. However, the converse is not necessarily true. For ex- 
ample, we find for the N-body models analysed below that 
the genus curve retains its random phase shape almost per- 
fectly at all smoothing scales, even when the smoothed den- 
sity field exhibits already strong non-Gaussian features that 
develop in the non-linear growth of structure. However, the 
amplitude of these genus curves is reduced compared to the 
expectation based on the power spectrum alone. 

This amplitude drop can be taken as evidence for phase 
correlations that develop in the weakly non-linear regime. In 
order to quantify this efi^ect, we define the amplitude drop 

(10) 



R = 



where Arp is the amplitude of the random-phase density field 
with the same power spectrum as the examined density field. 

There are difi^erent possible ways to measure A^rp. For 
a fully sampled N-body simulation with periodic boundary 
conditions, perhaps the simplest way is to Fourier transform 
the original density field, randomize the phases in Fourier 
space subject to the condition Sk* = (5-fc, therefore Gaus- 
sianizing the field, and then transform back to real space and 
measure the genus again. This procedure will be applied to 
the simulations, that we study in Section 5. 

For the case of PSCz, there are non-trivial phase cor- 
relations between the mask and the density field due to the 
odd shaped survey volume, which spoil this procedure. In- 
stead of constructing a representation of the corresponding 
random phase density field, the estimate of N-r-p needs to be 
based on a direct measurement of the second moment of the 
power spectrum. 



For this purpose, Vogeley et al. (199'^) estimated the 
power spectrum directly, but this is also complicated by the 
presence of the mask, and in addition it requires a subtle 
integration of the measured spectrum to derive Arp. 

As an alternative we propose to measure the variance 
(Tt (A) of the underlying density field as a function of smooth- 
ing length. This allows an estimate of (fc^) via the logarith- 
mic slope of rj^{\): 

From this the amplitude Arp can be obtained using equa- 
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tion (^. Our estimator for the variance u^{X) is discussed 
in Appendix ^ In order to reduce the cosmic variance in 
the estimate of {k^') we actually fit the measured vari- 
ances with an analytic form based on a generic CDM power 
spectrum, and compute the derivative from this fit. 



2.4 Error estimates 

2.4-1 Uncertainty due to finite sampling 

We assume that the observed galaxy distribution represents 
a Poisson realization of the underlying density field. Since 
we have only one sampling realization of the PSCz density 
field at our disposal we employ the bootstrap technique to 
derive estimates of the uncertainties due to sampling noise 
in the PSCz genus results. 

For this purpose we generate an ensemble of bootstrap 
resamplings of the PSCz galaxy catalogue by randomly re- 
distributing the galaxies onto the sites provided by the orig- 
inal catalogue. Note that certain positions will be left empty 
while others will be occupied by several galaxies. An esti- 
mate for the statistical uncertainty of a measured quantity 
due to finite sampling may then be obtained as the variance 
among the measurements for the bootstrap ensemble. 

In the context of genus statistics this method has been 
proven to be a reliable and robust - if somewhat conser- 
vative - estimator of statistical errors. Moore et al. (1992) 



compared bootstrap errors for the genus statistic with the 
statistical errors obtained for different resamplings of full 
N-body simulations and found that the bootstrap errors are 
10 per cent higher than the 'true' statistical uncertainty. 

Of course, in addition to the sampling noise the genus 
curves are also uncertain due to cosmic variance resulting 
from the finite survey volume accessible to PSCz. We be- 
lieve that PSCz is the best sample examined so far in this 
respect and that due to the large number of resolved struc- 
ture elements the cosmic variance is smaller than in previous 
studies. However, perhaps the only reliable way to quantify 
its influence is to examine PSCz-like mock surveys extracted 
from large-volume N-body simulations. This technique al- 
lows a careful statistical study of systematic effects a nd wi ll 
be applied in a forthcoming paper by Springel et al. (1997). 



2.4-2 Shot noise bias 

In addition to the variance introduced by finite sampling 
there is also a systematic effect on the genus curve due to 
shot noise. 

In Appendix ^ we calculate an estimate of this bias. To 
the extent that the shot noise does not spoil the Gaussian 
nature of the underlying field, the genus amplitude may be 
estimated by 



(27r) 



+ 



„2 

"shot 



-1 3/2 



(12) 



Here, UBhot is the variance introduced by shot noise on a 
given scale, and is the variance in the field free from shot 
noise on that same scale. 

When the power spectrum follows a power law with 
slope n on the relevant scales, the amplitude becomes 



N 



(27r)2A3 



n + 3 



^2 

"shot 



(?) 



-1 3/2 



(13) 



We tested this relation with Monte-Carlo experiments 
and found it to be fulfilled very well. Since in a realistic 
application we typically have —3 < n < 0, shot-noise tends 
to increase the measured genus amplitude. 

In order to quantify the contribution of shot noise to 
the genus amplitude we introduce the ratio 

Aobs 



A 



Nt 



(14) 



where Aobs is the observed genus amplitude under the infiu- 
ence of noise and At is the amplitude of the field free from 
shot noise. 

Since the variance of the density field due to shot noise 
o"shot varies with distance, we compute an estimate of j4shot 
by averaging relation (^2|) over the survey volume. 

The variances and a^i^ot {r) are estimated as outlined 
in Appendix We also need the logarithmic slope (fc^) A'^ 
defined in equation (^l[). In practice, we want to obtain an 
estimate for ylshot that is not strongly affected by fluctua- 
tions due to cosmic variance between the different adopted 
survey volumes. We therefore base the estimate of A^hot on 
values for derived from an analytic fit to the measured 
variances, and on slopes {k'^') A'^ obtained for this fit. 

With these quantities computed, the integrations re- 
quired to estimate A^hot can be done. Typically, we find 
^shot ~ 1.1 — 1.2 for PSCz, so that genus amplitudes are 
increased by up to 20 per cent. It is therefore crucial to cor- 
rect amplitudes for shot noise if one is to draw quantitative 
conclusions from the measured genus amplitude. For this 
reason, we calculate a corrected amplitude N* = N/Ashot 
in order to estimate the amplitude drop R, which is partic- 
ularly sensitive to the value of the amplitude. 



3 CONSTRUCTION OF DENSITY MAPS 

3.1 The PSCz redshift survey data 

The sky coverage of the PSCz redshift survey is 84.1 per 
cent, excluding only the zone of avoidance, here defined by 
an infrared background exceeding 25MJysr~^ at 100 /im, 
and a few unobserved or contaminated patches at higher 
latitude. The excluded regions are coded in an angular mask 
as shown in Figure ^ 

We convert the observed redshifts to the Local Group 
frame and use them directly to infer comoving distances 
without further corrections for peculiar velocit ies. As has 
been demonstrated with N-body experiments ( Protogeros 



Weinberg 1997) and analytical calculations (Matsubara 



199C), redshift space distortions hardly affect genus statis- 



tics. Therefore this should be a good approximation for the 
genus analysis. However, a few blue-shifted galaxies had to 
be discarded, resulting in an effective sample of 14505 galax- 
ies. We will assume an Einstein-de-Sitter model for the back- 
ground cosmology throughout. The results should not be 
sensitive to this choice because in a cosmological context 
the PSCz density field maps only the local Universe. 
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Figure 1. Sky coverage of the PSCz survey in an AitofT projec- 
tion. The dots represent the galaxies in the survey, the shaded re- 
gions are unobserved and comprise the angular mask. The galactic 
center lies in the center of the plot. 



3.2 Smoothing procedure 

Assuming a universal luminosity function, an unbiased es- 
timate of the galaxy density field p{r) can be obtained by 
weighting the discrete point distribution m(r) of the ob- 
served galaxies with the inverse of the selection function 
S{r): 



p{r) (X 



m(r) 



(15) 



Here the selection function S{z) = {m{r)) is defined as the 
mean expected comoving number density of sources at red- 
shift z corresponding to the comoving position r. We employ 
the fitting form 



S{z) 



(16) 



and determine its parameters ( see T able ^) with the methods 
outUned in Springel & White (1997). Note that the selection 
function includes a correction for the strong evolution seen 
in IRAS galaxies. 

In order to obtain an estimate of the density field 
smoothed on some scale A we convolve p(r) with a Gaussian 
filter of the form 



W{x) = 



7r3/2A3 



exp --^ 



(17) 



Note that here we use Gott et al.'s (198!:) definition of a 
Gaussian filter rather than the conventional form. 

However, due to the lack of galaxies in the regions of 
the angular mask, the density would be systematically un- 
derestimated at locations close to unobserved patches if the 
smoothing were just done by a straightforward use of the 
kernel of equation (p^. In order to avoid this problem we 
employ the ratio method proposed by Melott & Dominik 
( 1993 ) , who have shown in a systematic study that a smooth- 



ing according to 



p(r) 



/ W{r -r')p{r')dr' 
J W{r -r")M{r")dr"' 



(18) 



leads to the smallest loss or distortion of topological infor- 
mation compared to a number of alternative schemes that 
treat the mask differently. Here M{r) is a mask field defined 
to be equal to for r lying behind the angular mask and 
to be 1 otherwise. For this choice the denominator of equa- 
tion (^^ essentially renormalizes the smoothing kernel to 
the survey volume visible from the reference point r. 

In the actual computation of the genus curve we only 
use the volume with M{r) = 1 which is not hidden by the 
mask. Additionally we restrict the genus computation to a 
sphere carved out of the smoothed density field. Note that 
there is no boundary smoothing effect due to the outer sur- 
face of this sphere since we also include the sources outside 
this final region in the smoothing process. 

The method we apply here has the advantage that it 
does not require some form of filling of the unobserved re- 
gions. As Melott & Dominik ( 1993 ) have shown, simple 
forms of such fillings that come to mind, like a constant den- 
sity padding or randomly placed points of the mean back- 
ground density, lead to larger biases in the genus curve than 
the ratio method. 

We compute the convolutions that appear in the numer- 
ator and denominator of ( p^ ) using a Fast Fourier Transform 
(FFT) on a 128^ mesh. We choose a grid size of 6 = A/6, 
which e nsures that the genus is free from finite mesh size 
effects (Hamilton, Gott & Weinberg 1986). The final depth 
J2max of the density field that we use for the topological anal- 
ysis is always small enough to avoid wrap around effects due 
to the periodic FFT smoothing. 



3.3 Depth of maps 

Because the galaxy density of a flux limited sample declines 
quickly with distance, the uncertainty of the smoothed den- 
sity estimate grows rapidly with redshift. It is desirable, of 
course, to use a survey volume that is as large as possible 
in order to beat down statistical noise and cosmic variance. 
However, with respect to the genus statistics the sampling 
must be at least dense enough to make discreteness effects 
in the genus curve negligible. 



Accordi ng to a useful rule of thumb ( Weinberg, Gott & 
Melott 1987) discreteness effects are small if 



X> d = S' 



(19) 



where d is the mean interparticle separation and S is the 
selection function. Adopting this criterion we choose a max- 
imal radius i?max given by A = ^(iimax)"^ and use it to 
delimit the usable survey volume K. This choice ensures 
that at the far edge of the survey volume the sampling con- 
dition is just met, and in the remainder of the volume the 
sampling is denser. 

However, if A is of the order of d one risks artificially 
introducing a meatball bias by identifying lonely tracers as 
isolated clusters. Therefore, any observed bias of this kind 
should be treated with caution, since it might well be a 
discreteness effect. 

Additionally, the genus amplitude is biased due to shot- 
noise present in the reconstructed density field, as we have 
discussed in Section 2. This noise contributes to second 
and higher moments of the density field. Unfortunately, the 
genus depends in a non-linear way on the various moments 



© 0000 RAS, MNRAS 000, 000-000 



The Topology of the IRAS Point Source Catalogue Redshift Survey 7 



400 



200 



100 



Table 1. Parameters of the selection function of PSCz. 
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Figure 2. The number of resolution elements (solid, left hand 
scale) for the PSCz survey when the maximal survey volume is 
used. Also shown is the radius (dashed, right hand scale) of the 
usable survey volume. 



of the density field, making it extremely difiicult to derive an 
unbiased genus estimator which is not affected by it. Hence, 
we estimate the contribution of shot noise to the genus signal 
a posteriori. 



3.4 Resolution elements 

The notion of number of resolution elements provides a use- 
ful way to roughly compare the statistical power of genus 
measurements. Because the smoothing extends over an ef- 
fective volume Vsm = n^^^X^ the number of independent 
structures that can be present in a finite survey volume is 
limited. This number is of order 



iVr, 



Vsm 



37r3/2A3 ' 



(20) 



where uj is the solid angle covered by the survey. 

The number A^'ros indicates the power of a data set used 
f or top ological analysis. With the QDOT survey Moore et al. 
( 19921 ) reached a maximum of about A^res = 80 whereas the 
most p owerful genus results so far came from Vogeley et al.'s 
(1994) treatment of the CfA survey, where they reached 
Nics ~ 260 for their best subsample. In a recent st udy o f 
the 1.2-Jy redshift survey Protogeros & Weinberg (1997) 
achieved A^ros = 170. 

The PSCz redshift survey can provide still more reso- 
lution elements as is evidenced in Figure ^ and Table ^. It 
also provides a high number of resolution elements over an 
unprecedented wide range of smoothing scales. In particu- 
lar, there are more than 300 resolution elements in the range 
6/i~^Mpc < A < 20/i~^Mpc. This wide dynamic range to- 
gether with the large volume covered make it more powerful 
than all previously examined samples. 

For the genus analysis we have examined smoothing 
lengths between 5/i~^Mpc and 56^~^Mpc with approxi- 
mately logarithmic spacing. In particular, all the smoo thing 
lengths consi dered in the studies by Mo ore et al. ( 1992 ), Vo- 
geley et al. ( 1994 ) and Rhoads et al. ( 1994 ) are contained 
in this set. Table ^ lists some relevant parameters for the 
different cases. 



a 



0.991 



0.02534 



h0.068 
-0.073 



-0.00130 
-0.00116 



3.445 



p 

+0.173 
0.158 



7 



1.925^ 



(141.3 ±2.4) X lO"'' 



Table 2. The smoothing lengths adopted for the topological anal- 
ysis of the PSCz survey. Listed are the adopted survey depth 
fimax, the resulting number N^cb of resolution elements and the 
number A^g^i of galaxies inside the survey volume. 

A 



h-iMpc] 


Rmax [/l-^Mpc] 






5 


34.92 


215.5 


2295 


6 


47.16 


307.3 


3550 


7 


58.35 


366.3 


4928 


8 


68.58 


398.5 


5909 


10 


86.91 


415.3 


7510 


12 


103.19 


402.3 


8775 


14 


118.04 


379.2 


9681 


16 


131.81 


353.7 


10356 


20 


156.96 


305.8 


11309 


24 


179.79 


266.0 


11968 


28 


200.97 


233.9 


12496 


32 


221.79 


210.6 


12833 


40 


263.44 


180.7 


13339 


48 


305.09 


162.4 


13684 


50 


315.50 


158.9 


13748 


56 


346.73 


150.2 


13924 



4 PSCZ RESULTS 
4.1 PSCz genus curves 

Figure ^ shows the genus curves of psczQ for smoothing 
lengths of 6, 8, 12, 16, 24, 32, 48, and 56/i"^Mpc. The 
other smoothing scales examined show similar looking genus 
curves. In each panel the solid line shows a high resolution 
curve computed for spacing = 0.01. While it is clear that 
adjacent points on the genus curve are highly correlated, the 
amount of jitter present in the curves nicely indicates the 
level of noise contained in the measurements. Also shown 
as a dashed line is the best-fit random-phase genus curve, 
which we use to infer the measured genus amplitude. 

We have chosen not to show the average genus curve 
over the set of bootstrap realizations (which would be a 
much smoother curve), because of the high degree of cor- 
relation between the different points in the curve. Further- 
more, the bootstrap generated curves are expected to be 
biased towards a meatball topology. However, we use the 
bootstrap generated curves to estimate the statistical un- 
certainty of individual points on the genus curve. Similarly, 
we use the bootstrap curves to estimate uncertainties for the 
genus meta-statistics. For this purpose, we estimate the un- 
certainty of a particular meta-statistic as the rms-fluctuation 
of the values obtained for the set of 15 bootstrap resamplings 
of the data. 



* The PSCz 
genus curves may be retrieved electronically at http://www.mpa- 
garching.mpg.de/PUBLICATIONS/DATA/971121_pscz 
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Table 3. Summary of numerical results obtained for the genus meta-statistics of PSCz. Listed for each smoothing length are the genus 
amplitude in the dimensionless form {2'rT)^X^N, the absolute genus amplitude G = ^NR^^^ in the survey volume, the shift Ai/, and 
the width Wi,. We also give the measured variance O"^ of the smoothed density fields, and our estimated shot-noise correction factors 
^shot '^'nd the amplitude drop R. 







Q 




Au 








0-2 




^shot 




R 




5 


0.20 ±0.03 


7.40 ± 1.12 


-0 


24 ±0.12 


2 


24 ±0 


20 


0.6691 ±0 


1312 


1 


09 





58 ±0 


09 


6 


0.31 ±0.03 


16.05 ± 1.44 


-0 


14 ±0.07 


2 


01 ± 


10 


0.6044 ± 


0772 


\ 


10 





81 ±0 


07 


7 


0.36 ± 0.04 


22.05 ± 2.38 


-0 


04 ±0.07 


2 


13 ± 


12 


0.5163 ±0 


0485 




10 





87 ±0 


09 


8 


0.41 ±0.02 


27.11 ± 1.23 





05 ±0.05 


2 


19 ±0 


17 


0.3937 ±0 


0315 


1 


11 





92 ±0 


04 


10 


0.45 ± 0.03 


31.20 ± 1.99 





05 ± 0.03 


2 


16 ±0 


08 


0.2463 ± 


0170 


1 


13 





92 ±0 


06 


12 


0.47 ±0.04 


31.58 ± 2.52 





02 ±0.05 


2 


14 ±0 


09 


0.1873 ±0 


0201 


1 


14 





88 ±0 


07 


14 


0.48 ± 0.03 


30.43 ± 2.18 


-0 


05 ± 0.05 


2 


11 ± 


10 


0.1474 ±0 


0156 


1 


16 





84 ±0 


06 


16 


0.51 ±0.05 


30.44 ±3.19 


-0 


03 ± 0.04 


2 


28 ± 


10 


0.1162 ±0 


0122 


1 


17 





84 ±0 


09 


20 


0.64 ±0.05 


32.75 ± 2.73 





05 ± 0.06 


2 


08 ± 


14 


0.0751 ±0 


0079 


1 


19 





94 ±0 


08 


24 


0.83 ± 0.06 


36.89 ± 2.59 





04 ±0.05 


1 


86 ± 


12 


0.0482 ± 


0055 


1 


20 


1 


12 ±0 


08 


28 


0.83 ± 0.07 


32.73 ± 2.59 





07 ±0.07 


1 


94 ±0 


18 


0.0336 ± 


0042 


1 


20 


1 


06 ±0 


08 


32 


0.79 ± 0.09 


27.97 ±3.19 





03 ± 0.07 


1 


81 ± 


14 


0.0239 ± 


0033 


1 


20 





95 ±0 


11 


40 


0.71 ±0.09 


21.53 ± 2.68 


-0 


06 ± 0.06 


1 


95 ± 


22 


0.0145 ±0 


0023 


1 


18 





79 ±0 


10 


48 


0.61 ±0.08 


16.54 ± 2.16 


-0 


01 ±0.08 


2 


11 ± 


20 


0.0103 ±0 


0019 


1 


16 





64 ±0 


08 


50 


0.68 ±0.11 


18.22 ± 2.93 





03 ± 0.09 


2 


06 ±0 


20 


0.0100 ±0 


0019 


1 


15 





72 ±0 


12 


56 


0.66 ±0.11 


16.65 ± 2.83 





06 ± 0.09 


2 


09 ± 


18 


0.0089 ± 


0018 


1 


12 





68 ±0 


12 
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Figure 4. The shift measured for the PSCz genus curves. The 
error bars are based on bootstrap resa mplin gs of the data. Also 
shown are the results of Vogeley et al. (1994) for the CfA survey. 
Here we indicate the uncertainty of CfA by the errors the authors 
report for mock catalogues extracted from a LCDM model. 



Figure 5. The width of PSCz genus curves, with error estim ates 
derived with the bootstrap technique. Vogeley et al.'s (1994) re- 
sults for the CfA catalogue are shown for comparis on. H ere the 



errors are taken to be the uncertainty Vogeley et al. ( 1994 ) report 
for mock catalogues extracted from a LCDM model. 



As is evident from the plots in Figure ^, the PSCz genus 
curves follow the random-phase expectation rather well. In 
particular, no obvious deviations like shifts, broadenings or 
the like are observed. 

This is also confirmed by Figures ^ and ^, which plot 
the meta-statistics shift and width for scales ranging from 
5h~^Mpc to 56/i~^Mpc. The error bars give the rms- 
fluctuations obtained over 15 bootstrap resamplings of the 
data. Acomparison is also given with the results of Vogeley 
et al. ( |l99^ ) for the CfA survey. 

Apart from a slight indication of a meatball-shift and 



a broac ening of the genus curves at small scales, these two 



meta-statistics fail to show any significant departures from 
Gaussianity. 

We turn now to the amplitudes of the genus curves 
themselves. Figure |^ displays our results for PSCz, together 
with the amplitudes of variants of CDM models, which we 



will discuss fully in the next section. Also shown is the ex- 
pected amplitude in linear theory for the power spectrum 
with shape parameter F = 0.2. The amplitudes of 
the simulations are based on the dark matter distribution 
which we identify with the galaxy distribution. It should be 
stressed that the use of volume weighting makes this a rela- 
tively weak assumption since we only need regions of higher 
smoothed galaxy density to correspond to regions of higher 
smoothed mass density but not galaxy and mass densities to 
be proportional to one another. The use of volume weig hting 
makes the ge nus curves insensitive to details of bias ( Park 



Gott 19911 ). 

Based on Figure |^, the standard cold dark matter model 
in which the galaxies trace the mass (SCDM) is clearly ruled 
out at high significance level. This can be traced back to the 
fact that the shape of its power spectrum is in strong dis- 
agreement with PSCz. The other three N-body models do 
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Figure 6. Genus amplitudes for PSCz and for the four CDM 
models. The error bars for PSCz represent bootstrap estimates 
of the sampling uncertainty. The solid line gives the expectation 
based on linear theory for a power spectrum with shape parameter 
r = 0.2. While PSCz matches the tCDM, ACDM, and OCDM 
simulations well, the SCDM model exhibits genus amplitudes that 
are too large. 
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Figure 7. The variance a'^(A) found in Gaussian cells for the 
PSCz survey. The Icr uncertainties are estimated from a set of 
PSCz-like mock catalogues extracted from a N-body simulation. 
The solid line represents the best fit of the power spectrum ( |2l[ ) 
to the data, with T = 0.19 ± 0.04 and o-g = 0.84 ± 0.04 For 
comparison also shown is the result of Saunders et al. ( 1991 ) for 
the QDOT survey (with their error bars). 



much better in this respect; they match the PSCz ampli- 
tudes very well, because the shape of their power spectrum 
is in good agreement with PSCz. 

It is important to point out the significance of these re- 
sults. Although a counts-in-cell or power spectrum analysis 
of the PSCz redshift survey will probably lead to the same 
result for the shape of the power spectrum, the route we fol- 
lowed here is radically different. The genus depends on the 
spatial coherence of structures rather than on the rms ampli- 
tude of fluctuations. Moreover, it depends only on the rank 
order of density values, not on the values themselves. In the 
presence of unusual forms of biased galaxy formation, non- 
Gaussian initial conditions, or observational errors, there is 
no guarantee a priori that the determination of the shape 
of P(k) via the genus amplitude will be consistent with that 



Figure 8. The amplitude drop R measured fo r PSC z. Also shown 
are the results obtained by Vogeley et al. ( 1994 ) for the CfA 
survey . Her e the errors are taken to be the uncertainty Vogeley 
et al. (1994) report for mock catalogues extracted from a LCDM 
model. 



infered from a direct measurement of P{k). The fact that 
the same conclusion is reached by a very different analysis 
of the data lends considerable, largely independent support 
to the P{k) result. 

Note that, in Figure ^ we have plotted the amplitude 
A'^ of PSCz without a shot-noise correction, i.e. these values 
are likely to be biased slightly high. 



4.2 Variance and amplitude drop 

Figure ^ shows our results for the variance of the Gaus- 
sian smoothed PSCz density fields, obtained with the unbi- 
ased estimator described in Appendix ^ The error estimates 
are based on an ensemble of PSCz-like mock surveys drawn 
from a N-body simulation. Note that we find some what more 
power on s mall scales than was found for QDOT (Saunders 
et al. 199l| ). 

A fit of the linear CDM power spectrum (^l|) to the mea- 
sured variances results in a shape parameter F = 0.19 ±0.04 
and a normalization erg = 0.84 ± 0.04, here expressed in 
terms of the rms-fluctuations in top-hat spheres of radius 
8/i~^Mpc. Note that these values refer to redshift-space 
only; no correction for redshift space distortions has been 
attempted. The quoted uncertainties are estimated by ap- 
proximating the measurements for different A as being inde- 
pendent. 

As outlined in Section 2.3.4 we also compute estimates 
for the logarithmic slope ofo^(A) and the approximate 
shot- noise correction factor Ashot- Combining these mea- 
surements allows us to estimate the amplitude drop of PSCz 
asR = N*/N,p = A/(yl3hotA,p). 

The uncertainty of the measured values for R is difficult 
to estimate, since we deal with a ratio of correlated quanti- 
ties. We tentatively assign the relative error of the amplitude 
measurement as error of R. Our experiments with mock cat- 
alogues suggest that this gives approximately the right size 
of uncertainty. 

In Figure ^ we plot our results for the amplitude drop 
of PSCz. Although some amount of phase correlation seems 
to be present at small scales, we find that PSCz does not 
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Table 4. Parameters of the examined CDM models. The simula- 
tions have been done by the Virgo collaboration. 





SCDM 


rCDM 


ACDM 


OCDM 


Number of particles 


256^ 


256^ 


256^ 


200^ 


Box size[ /i^^Mpc] 


239.5 


239.5 


239.5 


239.5 


^start 


50 


50 


30 


119 


no 


1.0 


1.0 


0.3 


0.3 




0.0 


0.0 


0.7 


0.0 


Hubble constant h 


0.5 


0.5 


0.7 


0.7 


r 


0.5 


0.21 


0.21 


0.21 


erg 


0.60 


0.60 


0.90 


0.85 



exhibit strong phase correlations on scales above 10 /i~^Mpc, 
in contrast to the findings of Vogeley et al. ( 1994 ) for the CfA 
survey. This means that the PSCz density field is consistent 
with random-phase initial conditions. In the next section we 
will further analyse this amplitude drop and compare it to 
the results of N-body simulations. 



5 N-BODY SIMULATIONS OF CDM MODELS 
5.1 The simulations 

We study four variants of the cold dark matter cosmogony. 
These models were kindly provided by the Vir go collabo- 
ration (Jenkins et al. 1997; Colberg et al. 1997) who have 
recently started an ambitious project to model the forma- 
tion of structure in the Universe making use of the largest 
available supercomputers in Europe. The simulations have 
been performed with an AP^M-SPH code named HYDRA 
(Couchman, Thomas & Pearce 1995). 



All four models we examine contain only cold dark mat- 
ter in periodic boxes of size 239.5 h~^Mpc. The basic param- 
eters of the individual runs are given in Table ^ Here the 
shape parameter F refers to the linear theory power spec- 
trum 

Pik) = — ^, (21) 

where a = 6.4 F"^ /i~^Mpc, b = 3.0 F~^ /i~^Mpc, c = 



1.7F-^ /i^^Mpc and ly = 1.13 


Bond & Efstathiou 1984 


Efstathiou, Bond & White 1992 


. The models are normal- 



ized so as to match the observed abundance of rich clusters 



of galaxies (Eke, Cole fc Frenk 1996) 



The standard cold dark matter model (SCDM) has the 
power spectrum (|2^) with a shape parameter F = 0.5 in the 
linear regime, as predicted from a CDM inflationary sce- 
nario with a primordial scale-invariant Harrison-Zel'dovich 
spectrum. The rCDM model is a variant of CDM with more 
power on la rge scales, as suggested by numerous recent ob- 



servations (Saunders et al. 1991; Peacock & Dodds 1994 



Oliver et al. 1996). The same shape of the power spectrum 
may also be obtained in low-density universes. We consider 
two models of this class; one open universe (OCDM), and 
one with a cosmological constant (ACDM) that retains a flat 
background geometry. 

The simulations contain typically more than 16.7 mil- 
lion particles, thus providing excellent data which allow very 
accurate measurements of the genus statistic, especially on 
small scales, where the number of resolution elements is high 




10 15 

X I /I 'Mpc] 



Figure 10. The amplitude drop measured for the four N-body 
simulations as a function of the smoothing scale. Also shown are 
the values estimated for the PSCz survey (and error bars). 



and the genus curves are hardly influenced by cosmic vari- 
ance. This is, however, not true for scales > 20ft~^Mpc. 
The number of resolution elements in the simulations drops 
from 19861 for A = 5/i"^Mpc to 180 for A = 24/i-^Mpc. 
Hence the survey volume is quite limited even for such large 
simulations. In fact, beyond 20/i~^Mpc the volume probed 
by PSCz becomes larger than that of the simulations. We 
therefore restrict our genus analysis of the models to the ten 
smoothing scales between 5/i~^Mpc and 24/i~^Mpc used in 
the analysis of PSCz. 

In order to construct smoothed density flelds for the 
simulations we bin the particles on a 128^ grid by CIC as- 
signment and subsequently smooth with a Gaussian kernel. 
We also compute Gaussianized representations of the simula- 
tions by taking the Fourier transforms of the original density 
fields, randomizing the phases in Fourier space subject to the 
reality constraint 5k* = 5-k and transforming back to real 
space. These fields are then smoothed in the same way as 
the original ones, and the genus curves are calculated as be- 
fore. However, due to the limited simulation volume there is 
some variance introduced in the resulting genus curves from 
the variety of different possible random-phase realizations. 
In order to reduce this effect to a minimum we compute 
10 Gaussianized flelds and take the average genus curve of 
those as the Gaussianized genus curves. The amplitude of 
the latter is used to estimate the amplitude drop. 



5.2 Results 

Figure ^ shows representative genus curves obtained for the 
SCDM and ACDM models. In each panel the heavy solid 
line gives the genus curve for the density fleld of the evolved 
simulation at the indicated smoothing scale. A fit of these 
data to the generic random-phase genus curve is shown as 
a dashed line. The dotted line shows the genus curve for 
the corresponding Gaussianized field, and the thin solid line 
marks the genus of the initial conditions of the simulations. 
For the sake of brevity we don't show further genus curves 
for the rCDM and OCDM models. They look qualitatively 
very similar to those in Figure ^. 

Note that there is no sampling noise in these density 
fields; the residual tremble that is present in the N-body 
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Figure 9. Genus curves of the SCDM (left column) and ACDM (right column) simulations at selected smoothing lengths. In each panel 
we show the actual genus curve of the evolved simulation as heavy solid line, and we give the best-fit random-phase genus curve as 
dashed line. The dotted line is the genus of the corresponding Gaussianized field, while the thin solid line shows the genus curve of the 
initial conditions. 
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Figure 11 . Thr; amplitude drop revisited. Here we plot the data 
of Figure hu again, however this time as a function of the variance 
of the smoothed density fields. The good match of the three F = 
0.21 models in this representation suggests that the differences 
that show up in Figure between these models are only due 
to their slightly different normalizations. The SCDM model (F = 
0.5), however, continues to show a considerably smaller amplitude 
drop. 



genus curves is caused by the finite number of resolution 
elements, i.e. ultimately by cosmic variance. 

It is apparent that all the genus curves follow the Gaus- 
sian w-shape very well. Apart from the small bubble-shift of 
the A-model at the smallest smoothing scales the only meta- 
statistic able to measure a significant deviation from Gaus- 
sianity is the amplitude. It is interesting to note that the 
Gaussianized genus curves of the SCDM model remain close 
to the ones of the initial linear density field at all smoothing 
scales. In contrast, the ACDM model shows signs of non- 
linear evolution of the power spectrum at A = 5 ^~^Mpc. Its 
effect has been to increase the effective slope of the power 
spectrum (lenkins et al. 1997) on the relevant scales, and 
thereby also to increase the genus amplitude of the Gaus- 
sianized fields. Higher-order correlations, however, 'keep the 
genus in place', and decrease the amplitude at small A below 
the linear theory expectation. 

We defined the amplitude drop as the ratio of the genus 
amplitudes of the evolved density field and its Gaussianized 
counterpart. Hence, this quantity is a quantitative measure 
of higher-order correlations. In all of the four models an 
amplitude drop is present at small scales, being largest for 
the ACDM model and smallest for the SCDM model. 

In Figure |lo| we plot this amplitude drop against the 
smoothing scale for the four simulations. As expected, the 
amplitude drop becomes smaller with increasing smooth- 
ing length. This is consistent with the picture that all non- 
linear features are sufficiently smoothed out in this regime. 
It is interesting however, that the models exhibit a varying 
strength of the amplitude drop. For the ACDM and OCDM 
models, phase correlations up to the scale of 24/i~^Mpc are 
detected, whereas the SCDM becomes very close to Gaussian 
at about 10/i~^Mpc, and the rCDM takes an intermediate 
position. 

Figure |l^ shows that part of this difference is caused 
by the different normalizations of the models. Here we have 



the density fields. The good match between the models with 
shape parameter F = 0.21 suggests that the strength of the 
amplitude drop at fixed a is primarily governed by the shape 
of the power spectrum. 

In Figures |l^ and ^ we have also plotted the measured 
amplitude drops for the PSCz survey. It is encouraging that 
the measurements for PSCz exhibit about the right strength 
of phase correlations like the one expected in CDM scenar- 
ios. In principle, one could have hoped that this measure- 
ment might be able to strongly favour one of the CDM vari- 
ants. However, at this point the discriminative power of the 
genus test is not strong enough for this. The ACDM, rCDM, 
and OCDM fit the PSCz almost equally well. Only Figure 
[l^ suggests a slight preference for rCDM, essentially favour- 
ing its normalization. Given the considerable uncertainties 
in the biasing issue, this is not more than circumstantial ev- 
idence for this model. Note that in contrast to Figure ^ the 
horizontal positions of the PSCz points in Figure hll depend 
on bias; a positive bias factor would reduce the mismatch 
between PSCz and the models at large . 

It should be noted, however, that we are comparing 
observations with perfect simulation data. Mock catalogues 
should be used to accurately assess systematic effects from 
shot noise, redshift space distortions, selection function, and 
so forth, so that one can establish firm conclusions about 



the viability of models. Protogeros & Weinberg (1997) have 
pioneered this approach, and it will also be applied in a 
forthcoming paper by Springel et al. (1997). 



6 CONCLUSIONS 

We have analysed the topology of a new large redshift sur- 
vey, the PSCz, by means of genus statistics. The genus test 
provides not only information about the shape of the power 
spectrum but is also sensitive to higher order correlations 
of the density field. The PSCz survey is well suited for the 
genus test, due to its large survey volume and its near full 
sky coverage. In particular, it allows a topological analysis 
with larger statistical significance than any previously ex- 
amined sample. 

The genus curves of PSCz are featureless at smooth- 
ing scales ranging from 5/i~^Mpc to 56/i~^Mpc, i.e. they 
exhibit the w-shape that is characteristic of random-phase 
density fields. In particular, we find no clear evidence for 
shifts or broadenings of the genus curves. Evidence for sig- 
nificant non-Gaussian signatures on large scales found pre- 
viously in Vogeley et al. (1994) are not confirmed by PSCz. 



When the genus amplitude is examined, we find that 
the PSCz density field is consistent with a CDM power spec- 
trum with shape parameter F ~ 0.2 and inconsistent with 
the SCDM model. These results seem consistent with pre- 
vious anal ysis of the shape o f the power spectrum based on 
topology ([Moore et al. 199^ ; Vogeley et al. 1994 ; Protogeros 
& Weinberg 1997). We expect that an ordinary power spec- 



plotted the amplitude drop against the variance a (A) of 



trum analysis of PSCz will give a similar result. Again, we 
want to stress that these two methods of estimating the 
shape of the power spectrum are largely independent. 

We have demonstrated that the genus curves of CDM 
models retain the w-shape of the random phase genus curve 
even at small smoothing lengths where non-linear gravita- 
tional clustering has already generated significant skewness 
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of the one-point PDF. However, the non-hnear evolution 
manifests itself in a depressed genus amplitude compared 
to the expectation based on the power spectrum alone. This 
amplitude drop directly quantifies the amount of higher or- 
der correlations present in the density field. When plotted 
against the variance of the smoothed density field, we find 
that the amplitude drop depends on the shape of the power 
spectrum alone. 

The measured amplitude drop for PSCz is consistent 
with the CDM models, i.e. the non-detection of strong phase 
correlations in PSCz for smoothing scales above 10/i~^Mpc 
supports the hypothesis that structure grew from random- 
phase initial conditions. Alternative models for structure for- 
mation that provide a non-Gaussian seed field (for example 
by means of cosmic strings or explosions) are expected to 
exhibit stronger higher-order correlations in their density 
fields. Future analysis of the genus of such models should 
be used to test whether these correlations indeed show up 
as a strong amplitude drop or a distortion of the shape of 
the genus curve. In this case, our results for PSCz can place 
strong limits on the viability of such models. 
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APPENDIX A: CLUSTERING AMPLITUDE 
DETERMINATION 

In this appendix we outline our method to calculate the 
variance cr^{X) of the smoothed PSCz density field. This 
measurement is used to estimate the shot noise influence on 
the genus amplitude and to obtain an estimate of the genus 
amplitude of the Gaussianized version of the P SCz density 
field. Our approach is that of Springel & White (1997) which 



is a slightly mo dified version of the method proposed by 
Saunders et al. (|l99l|). 



We start with an unsmoothed density field derived by 
binning the galaxies on a fine mesh. We assume that the 
probability of finding a galaxy in a voxel i of volume SVi 
is given by p cx piSiSVi where pi is the density and Si the 
value of the selection function in the cell. For convenience, we 
assume p = 1 in the following. Then the probability density 
function that describes the distribution of the counts in cell 
i is given by 



Pi{m) 



Erui 

JV=0 



(Al) 



Here roi = piSi denotes the expected number of galaxies in 
the cell. Moments of the counts in this cell can be derived 
from the moment generating function 



Mmi{t) — i Pi{m) e™'*dm = exp [mi(e* — 1)] 



by differentiating 



dt" 
at zero lag. 

We now focus on the smoothed density field 



di 



(A2) 



(A3) 



(A4) 



where Wij is the value of the effective smoothing kernel (i.e. 
taking into account the renormalization due to the ratio 
method) between cells i and j. 

Because the rnT are independent variables with respect 
to the sampling process the moment generating function for 
di is simply 



MdAt) = Y[Mrr., [^t] =exp 



Hence, di is an unbiased estimate of the smoothed un- 
derlying field Pi = '^jWijpj- An unbiased estimate of the 
mean density d is therefore given by 



d = 



(A5) 



for arbitrary weights gi. We choose gi = [var(di)] ^ which 
provides a minimum variance estimator for d. To compute 



di' 



(A6) 



we first average over sampling realizations (overbar) and 
then with respect to the density field pi (angular brackets). 
We now want to determine the moments 



Rn 



- 1)") 



(A7) 



and, in particular, R2 — a^. Introducing the moment gener- 
ating function 



mj(e 



1) - dit 



the central moments are given by 



D 



(n) 



(di ~ diY 



d"M, 



di" 



(A8) 



(A9) 



Note that due to the rescaling we have di = pi. We can now 
relate the moments of the measured field with those of the 
underlying density field by using a binomial expansion 



{d, - 1)" 



fc=0 



k 



(d,-d,)'(p«-l)' 



(AlO) 



In particular, for n = 2 we find [di — 1)^ = D^^^ + [pi — 1)^. 

Then an unbiased estimate for the variance of the 
smoothed underlying density field is given by 



R2 



(d, - If 



D 



(2) 



(All) 



for arbitrary weights hi. Again, we want to estimate R2 with 
minimum variance by choosing the weights as 

1 



hi 



(A12) 



In order to co mpute these weights we make again use of 
equation (AlC): 



{di 



Dr'+4Di'Hp.~ir 

+ &Df\p.-lf + {p. 



(A13) 



Here the moments are given by equation (A9) as 



D 



(2) 



5^ 



Pi g2 ' 



D 



(4) 



pj- 



(A14) 



We now approximate the sums in these expressions by 
replacing the density pj by the averaged (smoothed) density 
Pi in the particular region. We then obtain 

D'f^^p^Y^^ Df^^PiY^, 

i?r'=3(p.F<^')%p.F/^\ (A15) 
where F/"' is given by 



y(") _ 



\ - Wi 



(A16) 



Of course, this is only strictly correct for a homoge- 
neous density field, but in the context of this analysis it can 
be expected to be a good approximation. Note that even 
if the weights hi are obtained only approximately, the final 
estimate of R2 remains unbiased. 

We can now take the ensemble average of equation 
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( A13) and compute the weights hi. Here a well known prob- 



lem in counts-in-cells analysis arises: the determination of 
second and higher moments involves the whole hierarchy of 
moments. A minimum variance estimate of R2 requires the 
knowledge of and R4, which in turn depend on still higher 
moments. To close this hierarchy we assume that the den- 
sity field is close to Gaussian, which implies i?3 = and 
Ri = S-RI • With this assumption the weights become 



= 2 y: 



K2) 



+ 2Ri 



+R2 < 3 



(2) 



(2) 



(A17) 



An estima te of R2 can now be obtained by iteratively solving 
equation (All). 



Since the di are strongly correlated an estimate of the 
statistical error in the determination of R2 by directly com- 
puting the variance of R2 in the above formalism is unre- 
alistic. Instead, we estimate the uncertainty in R2 by using 
an ensemble of PSCz-like mock catalogues extracted from a 
N-body simulation. 



APPENDIX B: GENUS AMPLITUDE AND 
SHOT NOISE 

The genus amplitude of a Gaussian random field depends 
only on the behaviour of the two-point correlation function 



of the 



Wt! Lhtii'tiforti 



sjiuooLlitid dyuyity fltild iieai' Llit; origin. VVt! 
estimate {() ) uudei Llm iuflueuce uf shut uuise, alluwiug the 
noise level to vary across the survey volume. 

We start by considering a discrete field m(r) = 
y^ ^. (5'°-'(r — ri) of points arising in a Poisson sampling pro- 
cess of some underlying density field St{r). 

Allowing for a radial variation of the expected number 
density S{r) = (m(r)) of tracers, we define the measured 
density contrast as 

Sir) = ^ - 1. (Bl) 
S{r) 

This gives rise to a smoothed density field 

5{r) = j W{r -r')5{r')Ar (B2) 

with an autocorrelation function 

|(r) = (5(ro)5(ro+r)). (B3) 

Performing t he en semble average in equation (B3) with 
Feldman et al.'s (1994) relation 



{m{r')m{r")) = S{r')S{r") [l + ^.{r' - r")] 
+5(r')5<°^(r'-r"), 

we find 



(B4) 



W(r')W{r - 



r')^t{r") dr' dr" + 



W{ro — r')W{ro + r — r') 



S(r') 



dr', 



(B5) 



where ^t(r) is the autocorrelation function of the underlying 
field. 

The second term in equation (|B5|) describes the shot 



noise contribution which varies across the survey volume. 
Hence the dependence on ro remains even after the ensemble 
average. Assuming that S{r) varies sufficiently slowly over a 
smoothing length and choosing a Gaussian smoothing kernel 
for W{r) we can substitute S{r') by S{ro) in the integrand 
because W{ro — r') is strongly peaked at ro. We then obtain 



^(r)=Ct(r-)+?shot(r) 
where 

1 



(B6) 



(27r)3/2A3 



ft(r") exp 



{r-r"f 



2A2 



dr ' 



(B7) 



is the contribution of the underlying density field and 



Cshot(r') = 



exp 



r 

2A2 



(27r)3/2A35'(ro 
stands for the shot noise contribution. 



Equation (B7) reads in Fourier space like 
Pt(fc) = Pt(fe) exp ' 



(B8) 



(B9) 



Adopting isotropy for the underlying field, we can also write 

(BIO) 



kr 

We now assume that the density field can approximately 
be taken to be Gaussian. In this case the a mplitude of the 
genus curve depen ds only on ^(0) and ^"(0) (Hamilton, Gott 



& Weinberg 1986). The previous equations give the corre- 



sponding contributions from the underlying field as 



and 



Pt(fc)d^fc 



C"(0)= j -Pt{k)d'k = 

while the shot noise gives rise to 
1 



:eshot(0) = 



(27r)3/2A35'(ro) 



and 

Cshot(O) = 



1 



2 

""shot 



(BU) 



(B12) 



(B13) 



(B14) 



(27r)3/2A5S(ro) A^ ' 

We can now estimate the amplitude of the genus curve 
under the influence of shot noise. It is given by 

3/2 

(B15) 



N 



1 



(2^)2 
1 

(2^ 



-r(o) 



+ 



(k^y 



3/2 
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